DTP/97/32 
hep-ph/9705229 
May 1997 



iS' Hadronic Antenna Patterns as a Probe of 

Leptoquark Production at HERA 



in 



> : b 

o\ 

(N 
(N 

in 
o 

On 



Oh: 
a: 



X 



M. Heyssler a and W. J. Stirling a ' b 

a Department of Physics, University of Durham, Durham, DH1 3LE 
Department of Mathematical Sciences, University of Durham, Durham, DH1 3LE 



Abstract 

Hadronic antenna patterns can provide a valuable diagnostic tool for probing the 
origin of the apparent excess of high x, Q 2 events at HERA. We present quantitative 
predictions for the distributions of soft particles and jets in standard deep inelastic 
scattering eq — > eq events and in events corresponding to the production of a narrow 
colour-triplet scalar ('leptoquark') resonance. The corresponding distributions of soft 
photon radiation, which are sensitive to the leptoquark electric charge, are also pre- 
sented. The distribution for one particular leptoquark assignment is shown to contain 
a radiation zero. 



1 Introduction 



The observation of an apparent excess of deep inelastic scattering events in positron-proton 
collisions at high Q 2 by both the HI [Q and ZEUS @ collaborations at HERA has prompted 
much speculation about possible new physics explanations. Obvious candidates are a new 
four-fermion contact interaction A~ 2 eeqq with A ~ (9(1 — 2 TeV), or the production of 
a new heavy 'leptoquark' resonance e + q — > LQ — > e + q with Mlq ~ 200 GeV ||. The 
electric charge of such an object is not yet known, but if clq = +2/3, corresponding to 
e + d — > LQ for example, then the new particle could be a heavy squark in an i?-parity 
violating supersymmetric extension of the Standard Model. 

It is important to investigate all possible ways in which one could distinguish between 
a conventional explanation (i.e. a fluctuation of the standard model DIS process) and new 
physics scenarios. One diagnostic tool which has already been advocated [|J in the context 
of large Et jet production at hadron colliders is to use the pattern of hadronic energy flow 
in the event [0-0] • This is based on the idea || |6| that the overall structure of particle 
angular distributions in a hard scattering process (the 'event portrait' or 'antenna pattern') 
is governed by the underlying colour dynamics at short distance. Thus the event portrait 
can be regarded as a 'partonometer' |§ mapping the basic scattering process. 

The events at high x and Q 2 at HERA seem ideally suited to such a study, being char- 
acterized by an energetic, well-separated lepton and jet in the final state. Furthermore the 
various candidate underlying eq — > eq processes (t-channel 7*, Z exchange, contact interac- 
tion, s-channel colour-triplet resonance production, . . . ) have distinctive antenna patterns, 
as we shall see. In practice one could, with sufficiently high statistics, use an additional 
soft (gluon) jet as a probe of the antenna pattern. With fewer events the distribution of 
soft hadrons can be used instead. Both of these quantities are related to the inclusive soft 
gluon distribution in the next-to-leading order eq — > eqg processes, the former directly and 
the latter through the hypothesis of Local Parton Hadron Duality (LPHD) [TJ| in which 
the angular distribution of soft particles emitted at wide angles to the energetic jets follows 
that of the underlying soft partons, with the rate being determined by overall multiplicative 
energy-dependent cascading factors || [15 . 



The idea, then, is to use the angular distributions of soft particles or jets as a probe 
of new physics contributions to high-Q 2 e + q scattering. We imagine a situation where a 
larger sample of (presumably standard model DIS) events at slightly lower Q 2 is used as a 
control, to check the approximate validity of our quantitative predictions for the antenna 
pattern. This can then be compared with the observed antenna pattern for the sample of 
excess events. As we shall see, in some cases the 'signal' and 'background' distributions can 
differ by factors of 2 or more. The variation of the patterns with the DIS variable y will also 
be a useful discriminant. 

We should also remark that experimental support for the feasibility of such antenna 
pattern studies has come recently from the CDF and DO collaborations at the Tevatron pp 
collider |T(J I7| . Their analyses have shown that distinctive colour interference effects in 
multijet and PU+jet production survive the hadronization phase and are clearly visible in 
the data. 

In the following section we derive the basic antenna pattern results for standard DIS and 
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leptoquark production. The case of a new contact interaction is obtained as a limiting case 
of the latter. In Section 3 we present numerical predictions for various typical kinematic 
configurations. Analogous results for soft photon production are obtained in Section 4, and 
Section 5 contains our conclusions. 



Antenna patterns for deep inelastic scattering and 
leptoquark production 



The distribution of soft gluon radiation is controlled by the basic antenna pattern (see for 
example Ref. [0) 



Pi ■ Pj 1 - n* • rij 



Pi ■ k pj ■ k u 2 (l — n • rij) (1 — n • rij) ' 

where the pf = Ei(l, rij) are the four-momenta of the energetic quarks and leptons partic- 
ipating in the hard scattering process, and k^ = u(l,n) is the four-momentum of the soft 
gluon. The radiation patterns presented below correspond to the uj/Ei — > limits of the 
exact eq — *> eqg matrix elements. 

We start by considering the Standard Model process e(pi) + q{p-z) — > e(ps) + q{p^) by t- 
channel 7*, Z exchange. If the invariant mass of the eq system is M, and if the angle between 
the incoming and outgoing quarks (in the eq c.m.s. frame) is Q q , i.e. cosO g = n 2 • n 4 , then 
the usual DIS variables are 

M 2 1 

x = — , j/ = -(l-cose,), Q 2 =yM 2 . (2) 

S 

The scattering process with the various momenta labelled is shown in Fig. 1. 

Since our aim is to distinguish the patterns for resonance production and the normal 
deep inelastic scattering, we consider fixed M and variable Q q (y). For the Standard Model 
process the gluon energy and angular distribution is simply (see for example [0) 

1 da Ci s Cp 



a duidVL n 4n 2 



u -^SM (3) 



where 



T = 2124] = 2P2 ' P4 = 2(1 - cos 6 9 ) 

SM 1 J p 2 -kp A -k u 2 {l-cos6 2 ) (1-COS0 4 ) ' 
where cos 9{ = n ■ denotes the angle between the soft gluon and the corresponding quark. 
The gluon emission is coherent, and depends on the relative orientation of the incoming and 
outgoing quark directions. Eq. (ffl) can be interpreted as a colour 'string' connecting the 
incoming and outgoing quarks JT5|, |18j , and is closely related to the familiar result T = 2 [qq] 
for the crossed process e + e~ — > qq (see for example J7J). 

We now turn to the radiation pattern corresponding to the production of an unstable 
colour-triplet, s-channel scalar resonance LQ of mass M and decay width T, i.e. eq — > 
LQ — > eq. We first note that the emission of a soft gluon off an on-shell colour-triplet scalar 
boson is described by the same factor as emission off a colour-triplet fermion, i.e. 

* ^g.£j$& , (5) 
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where T a is a SU(3) colour matrix, P M is the momentum of the emitting particle, and e" 
is the gluon polarization vector. We can therefore use results already obtained for heavy 
quark production and decay P,|7|,|l9 H23ll to write down the result for leptoquark production 



and decay:[] 

Jl Q = 2 {[2P] + [4P] - [PP]) + 2 X ([PP] + [24] - [2P] - [4P]) (6) 
where P = pi + P2 is the leptoquark momentum. The factor x i n © is given by 

_ M 2 r 2 r 2 
x ~ (p ■ k) 2 + M 2 r 2 ~ cj 2 + t 2 ^' 

where the second expression corresponds to the LQ c.m.s. frame. As discussed at length 



in Ref. [pi]] , the radiation pattern depends, through the factor x-> on the relative size of the 
gluon energy and the leptoquark decay width. In this respect it is instructive to consider 
the two (formal) limits r —>■ oo (x — > 1) and r — > (x — *■ 0), for fixed u. In the former, 
the leptoquark decays immediately after it is produced and has no time to radiate gluons of 
wavelength ~ In this limit 

:Flq-2[24], (8) 

which is identical to the standard DIS pattern (f|) corresponding to coherent emission. In 
contrast, for T/u — > the emission takes place on two very different timescales, correspond- 



ing to the production stage and the decay stage |21[] : 



Jl Q -> {2[2P] - [PP]} + {2[4P] - [PP]} . (9) 

At threshold, where there is essentially no radiation from the heavy leptoquark, the two 
terms in {} correspond to independent radiation off the initial and final state (massless) 
quarks, see (|11]) below. Note that it is straightforward to verify that the first term on the 
right-hand side of (||) does indeed correspond to the k^ 1 — > limit of the real gluon emission 
matrix element squared for e + q — > LQ + g calculated in Refs. [^3], P5| . 

With no a priori knowledge of the decay width of the new heavy particle, the antenna 
pattern (Q) could in principle be used to obtain a measurement. This was the approach 
advocated in Ref. |21| for the top quark. As we shall see in the following section, in certain 
regions of phase space the antenna pattern is very sensitive to Xi an d therefore to V . In 
practice, it seems that for the class of leptoquark models proposed to explain the excess of 
high-Q 2 events at HERA, the decay width is rather small. In particular, a scalar leptoquark 
coupling with strength A to eq has a corresponding decay width T = M\ 2 /(16tc). For 'first 
generation' leptoquarks values of A 2 < O(10~ 2 ) are allowed by low-energy data (see for 
example Ref. |23j] and references therein). This implies that such resonances should be very 
narrow, i.e. T < O(40 MeV) for M ~ 200 GeV. If we are interested in the distributions of 
soft hadrons or jets with energies of order a few GeV, then x <S 1 and (^) is the appropriate 
distribution for the leptoquark signal. 

Finally, we note that the antenna pattern for a eeqq contact interaction corresponds to 
the limit x — > 1, and is therefore identical to the standard DIS result, Eq. (|4|). 

1 In fact the soft gluon distribution for eq —> LQ — ► eq is identical to that for Wb — » t — > Wb with 
m t = Af LQ , T t = T L q and m b = 0. 
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3 Numerical results 



In this section we present numerical results for the standard model DIS and leptoquark soft 
gluon distributions. We work in the eq c.m.s. frame with angles defined as in Fig. 1, and 
focus on the dependence of the dimensionless quantity Af = uj 2 ^, where JF SM and ^lq are 
defined in (|) and (|9]) respectively, on the gluon direction n. Simple algebra gives 

Af 2(1- cos 6 g ) 

(1 — cos u 2 ) {1 — COS #4) 



. 1 + cos 6 2 1 + COS 64 

A/lq = -r- + -r- , 

1 — cos 62 1 — cos 64 



Af L Q _ 1 - COS 62 cos 64 

A/cjM 1 COS Qq 

The patterns and their ratio are displayed in Figs. 2 and 3, as functions of 9 g and <ft g , the 
polar and azimuthal gluon angles with respect to the incoming quark direction,^ and for 
fixed values of Q q = 45°, 90°, 135°, i.e. y = 0.146,0.5,0.854. To avoid the collinear-singular 
regions of phase space, cuts 9 2 , 64 > 10° are imposed.^ 
We note the following points: 

(i) For the standard model distribution, there is a significant enhancement of radiation in 
the region between the quark directions (i.e. ~ 0°, 0° < 6 < Q q ), as expected. This 
enhancement is largely absent in the LQ case, where the radiation pattern is simply a 
superposition of independent radiation off the initial and final state quarks. 

(ii) In the limit Q q — > 0°, A/sm vanishes everywhere since the final state comoving colour 
triplet and antitriplet behave as a colour singlet, whereas A/lq is simply twice the 



radiation off a single quark. In Ref. [22|, similar effects where discussed for e + e 



tt — > W + W hh production at threshold. 

(iii) For Q q = 90° scattering, the ratio of the SM and LQ distributions achieves its minimum 
and maximum values in the plane of the scattering, thus A/"lq = |A/sm for 6 g ) = 
(0°,45°) and (180°, 135°) and A/" L q = |A/" SM for {<j> g , 6 g ) = (0°, 135°) and (180°, 45°). 
The distributions are the same for gluon directions in the planes perpendicular to n 2 
and n 4 , i.e. 6> 2 , O4 = 90°. 

Finally, from the above discussion we would expect that the azimuthal distribution of 
soft gluons (hadrons) around the final state quark (jet) direction would be more uniform 
for quarks from leptoquark decay than from standard deep inelastic scattering. To see this, 
we show in Fig. 4 the azimuthal 4> g distribution of the gluon around the final state quark 
direction n 4 , for Q q = 90° and various fixed 64. A significant azimuthal asymmetry for Asm 
is observed with a maximum in the plane of the scattering between the quark directions 
(4> g = 0°), as expected. In contrast, the dependence of Alq on cj) g is very weak, particularly 
for small 64. 



9 g sin 8 g + cos 9 g cos 6 9 . 



3 The cuts on 02, O4 are omitted in Fig. 3, since the ratios are finite (= 1) in the two collinear limits. 
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4 Soft photon emission 



As discussed in the introduction, it would be of considerable interest in distinguishing new- 
physics models of the HERA high-Q 2 events to know the electric charge of the quarks in 
the eq — > eq process. In principle, this information is contained in the distribution of soft 
photon radiation, which can be obtained in an analogous way to the soft gluon distributions 
of Section 2. The main difference is the presence of additional contributions from emission 
off the incoming and outgoing positrons]^ The result is (cf. Eqs. (P,§J|)) 

1 da a . 

cu 7 (13) 



ctq doj~dVL n An 2 



where 



^s 7 m = e 2 [24] - e q {[12] + [34] - [14] - [23]} + [13] (14) 
= e 9 (l + e J{[2P] + [4P]} + (l + e? ){[lP] + [3P]} 
-e q {[12} + [34}}-(l + e q ) 2 [PP] 
+ X ((1 + e q ) 2 [PP] - e q (l + e q ){[2P] + [4P]} 

-(1 + e q ){[\P] + [3P]} + e 2 [24] + e,{[14] + [23]} + [13]) (15) 

and, as before, J 7 g M = ^lq(x — !)• As argued in the previous section, it is the x "~ > limit 
of ^"lq which is relevant in practice, i.e. for photons with energy uj^ ^> I\q. In this limit 
we have 

= e g (l + e 9 ){[2P] + [4P]} + (l + e ,){[lP] + [3P]} 
-e ? {[12] + [34]}-(l + e 9 ) 2 [PP] 

= H(cos6 2 ) + P(cosfl 4 ) (16) 



where 



P^ = TT; + r ^-|(l + e 9 ) 2 . (17) 

An interesting feature of the above distributions is the presence of radiation zeros (see for 
example |2H||), i.e. directions of the photon three-momentum n for which the cross section 
vanishes. To see this for the distribution (^|) we note that 

H = for z = z = 1^1 . (18) 

1 + e q 

For the two cases of interest e q = |, — | for which z = |, 2. Therefore only for e + u scattering 
is the radiation zero in the physical region. For the full distribution ([II]) to vanish we 
obviously require 

cos 6*2 = cos 6*4 = Zq . (19) 

The results in this section are for e + q — > e + q scattering. Those for e~q — > e~q can be obtained by an 
appropriate change of sign. 
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Thus for e + u scattering the radiation zero is in the direction given by the intersection of the 
two cones of half-angle 0q = cos _1 (l/5) = 78.46° centred on the quark directions ri2 and 114. 
Three cases can be distinguished: 

(i) For 0° < Q q < 20 o there are two solutions, corresponding to 

«, = *, ^ ±C os-p«). (20) 

(ii) For Q g = 20 o there is one solution, 

0^ = O , 7 = 0° , (21) 
corresponding to the bisector of the quark directions in the scattering plane. 

(iii) For Q q = 0° there is a cone of solutions corresponding to # 7 = 0q. 

Although the above results on the location of the radiation zeroes have been derived for 
the leptoquark radiation pattern, they apply equally well for the standard model distribution 
(|H|), or indeed for the generic distribution flI5|) for arbitrary \- This follows from the fact 
that the zeroes are the result of completely destructive interference between the classical 
electric fields associated with the different charged particles, see for example Ref. JJ. They 
depend only on the relative orientation of the various particles, irrespective of whether 
intermediate resonances are formed. The key point to note is that since presumably the 
leptoquark couples to either u or d but not both, the radiation zero is either fully present 
or completely absent. In contrast, the standard model background is a linear combination 
(determined by the parton distributions) of u- and d-type distributions, giving a dip rather 
than a zero. 

As a numerical illustration of these results, we show in Fig. 5 the antenna patterns M^ui 
N2q{e q = 2/3) and J\f2c>( e q = ~ 1/3), with Q q = 90°. To exhibit the radiation zeroes more 
clearly, Fig. 6 shows the 7 dependence of the leptoquark e q = 2/3,-1/3 distributions at the 
critical polar angle # 7 = 0q, i.e. the slices through the two-dimensional distributions of Fig. 5 
at this value of # 7 . The two zeroes of the e + u distribution at the 7 angles given by Eq. ( |2"0D 
are clearly visible. Note also that the behaviour of the distributions near the positron and 
the quark jet directions simply reflects the magnitude of the charge of the corresponding 
particles. 



5 Conclusions 

If the observation of an excess of high-Q 2 events at HERA persists, it will be important to 
devise new analysis techniques for identifying the origin of the excess. In this paper we have 
shown that the angular distribution of the accompanying hadronic radiation - the antenna 
pattern - is a potentially powerful tool for discriminating standard deep inelastic scattering 
events from those arising from the production of a long-lived coloured scalar 'leptoquark' 
resonance. The main qualitative difference is the absence for the latter of an enhancement 
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of hadronic radiation between the incoming and outgoing quark jet directions (string effect), 
as shown in Fig. 2. It follows that soft hadrons are distributed more uniformly in azimuth 
around the final state quark jet direction in events where a leptoquark is produced, see Fig. 3. 
Our quantitative predictions are based on the phenomenologically successful principle of 
Local Parton Hadron Duality, and should therefore be a good guide to the behaviour of the 
distributions of soft hadrons and jets in the detectors. Ultimately, however, there will be no 
substitute for detailed Monte Carlo studies based on parton-shower/hadronization models, 
provided that these include the correct underlying colour structure. 

Finally we have extended our results to include soft photon radiation. Here the distribu- 
tions have an additional sensitivity to the electric charge of the leptoquark, which is a crucial 
parameter in distinguishing models. For the case of charge 5/3 leptoquarks, produced for 
example in e + u collisions, the soft photon distribution contains radiation zeroes. These are 
absent for charge 2/3 leptoquarks produced in e + d collisions. 
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Figure 1: Parametrization of the kinematics for e 4 ~(pi)q(p2) —> e+ G°3)<?G°4) + scattering 
in the e + q c.m.s. frame. The orientation of the soft gluon relative to the scattering plane 
is denoted by 9 g and (f> g or, alternatively, by the angles with respect to the directions of the 
participating quarks: # 4 and 9 2 = 6 g . 
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Figure 2: The dimensionless antenna patterns Msm = ^ 2 -7~sm [( a )>( c )>( e )] an d Alq = u 2 ^^ 
[(b),(d),(f)] of Eqs. ( pUjJTTI ) for different c.m.s. scattering angles Q q (cf. Fig. |TJ). Note the 
cut of 10° imposed around the incoming and outgoing quark directions. 
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Figure 3: The ratios A/lq/A/sm of the distributions in Fig. g for the three different c.m.s. 
scattering angles. In this case no angular cuts have been imposed. 
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Figure 4: The dependence of the antenna patterns Asm and Alq on the azimuthal angle 
<f)g of the soft gluon around the outgoing quark q{p^). The gluon direction describes a cone 
around the quark of half-angle 64. The direction of the incoming quark q{p2) is defined by 
(f)g = 0°, and the incoming positron e + (pi) is at <p g = ±180°. The overall c.m.s. scattering 
angle is fixed at <d q = 90°. 
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Figure 5: The pattern of soft 7 radiation according to Eqs. ( |i~4] , |T6|) with N"sm,lq = w r^sM,LQ' 
for e + d scattering [(a),(c)] and e + u scattering [(b), (d)]. The overall c.m.s. scattering angle 



is fixed at Q q = 90°. 
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Figure 6: The soft photon antenna pattern Alq for Q q = 90° at the critical angle 9 1 = 
9 ~ 78.46° for e + u (solid line) and e + d (dashed line) scattering. Note the radiation zeros 
at 7 ~ 78.22° (cf. Eq. (|20|)). The positions of the e + and q jets are indicated. Note also 
that J\f$ M shows quantitatively the same behaviour for this choice of 9 T 
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